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Abstract 

A T-shaped field-effect transistor, made out of a pair of two-dimensional 

electron gases, is modeled and studied. A simple numerical model is developed 

to study the electron distribution vs. applied gate voltage for different gate 

lengths. The model is then improved to account for depletion and the width 

of the two-dimensional electron gases. The results are then compared to the 

experimental ones and to some approximate analytical calculations and are 

found to be in good agreement with them. 
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I. INTRODUCTION 



Extensive work has been published on various nano-structure devices that use two- 
dimensional electron gases (2DEG)§1. Numerous devices have been built using nanofab- 
rication that involve shaping the 2DEG by applying gate voltages. Some of the techniques 
use in-plane gatingiB; however, most of them use off-plane gatingiHll. Here, we study an 
off-plane device whose gate is also a 2D EG. 

Our work was prompted by a paper by H. L. Stormer et al.@ In that paper, a new type 
of T-shaped field-effect transistor (T-FET) was presented with an ultra short gate length. 
The transistor is made out of a pair of 2DEGs that are oriented perpendicular to each other. 
One is used as the source-drain channel, and the other is used as the gate. When the gate 
voltage exceeds a certain critical value, a depletion region is created in the source-drain 
channel which cuts off the current flow through it. 

The device is made exclusively using molecular beam epitaxy (MBE). The 2DEG of the 

gate is made by standard modulation doping of a GaAs substrate in the following way: 
4.3 /im undoped Alo.3Gao.7As, 

200 A undoped GaAs quantum well, 

150 A undoped Alo.3Gao.7As spacer, 

1.5 x 10 16 m~ 2 Si 5 doping, 

20 jim undoped Alo.3Gao.7As. 

Following this, the structure is now cleaved on one edge, and the growth continues in a 

perpendicular direction. The growth sequence is as follows: 
200 A undoped Alo.3Gao.7As separator, 

150 A undoped GaAs quantum well, 

150 A undoped Alo.3Gao.7As spacer, 

2 x 10 16 mr 2 Si 5 doping, 

1700 A undoped Alo.3Gao.7As, 

600 A undoped GaAs. 

This produces a modulation-doped quantum well which is 200 A wide for the gate and 150 A 
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wide for the source-drain channel. A cross section of the transistor is shown in Fig. [|. 

In Sec. m we develop a very simple model for the transistor, based on work done by X. 
Liu and one of uM. In this model, we assume that the electron distribution is related to the 
potential energy by the Thomas-Fermi equation 

n = (5-U)u. (1) 

Here, n is the extra electron distribution due to the applied gate voltage, v = is the two- 
dimensional density of states, and 5 is the change in chemical potential of the two 2DEGs 
due to the bias. The equation is valid everywhere except inside the insulating barriers and 
any created depletion regions. It does not give an exact solution to the problem, but it 
can give a good understanding of the general behavior of this kind of semiconductor device. 
Furthermore, we initially assume that the two 2DEGs have no width in the direction in 
which they are confined by the quantum wells, and that the donors that created the gases 
are in the same position as the 2DEGs. We also assume that no depletion region is created. 
These assumptions allow us to come up with an analytic solution to the problem using 
Fourier transforms. 



In Sec. UJ, we allow for a depletion region to be created. The solution using Fourier trans- 
forms is no longer feasible, and a numerical solution is used throughout. Special techniques 
had to be developed for this. 

In Sec. |V[ we make a more realistic model of the transistor by taking into account the 
finite width of the gases. The techniques used to solve these are very similar to those used 
in Sec. 

In Sec. [V], an approximate analytic approach is used for the problem inspired by Ref. |9|. 
The results of the two methods are then compared in the conclusion. 



II. THE ONSET OF DEPLETION 



As a first approach, we assume that the two 2DEGs have no width and that no depletion 
region is created. Furthermore, we assume that the donors are in the same position as the 
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2DEGs. As shown in Fig. |2|, we choose the x axis for the gate and the y axis for the source- 
drain channel. The above assumptions simplify the equations for the potential energy of the 
gases to 



U(x,y) 
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which on each axis gives 
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In addition, for each of the 2DEGs the Thomas-Fermi equations give 

n(x) = [5 X - U(x)]u, 



and 



n(y) = [S y - U(y)}u 



(2) 



(3a) 



(3b) 
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for the x and y axis respectively. Substituting in for the potential energy terms on each axis 
from Eqs. (Q), and rescaling a few of the parameters we get 
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The new units used are: length: ao = ; density: 



(5b) 
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are two integral equations with two unknown functions n x (x) and n y (y). In addition to the 



above equations, the two additional requirements are the neutrality 

/*oo /*oo 

/ dxn x (x) + / dyn y (y) = 0, 

iff J — oo 

and the fact that 
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The above equation states the fact that the difference in the chemical potentials of the two 
sides is just the gate voltage. 

We can obtain one equation for n x (x) by using Fourier transforms to eliminate n y (y). 
The resulting integral equation for n x (x) is 
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In a very similar way, we can obtain an equation for n y (y) in terms of n x (x). By manipulating 
Eqs. (|5|), we can express the Fourier transform of n y (y) in terms of n x (x) as n y (k y ) = 
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which, after taking the inverse Fourier transform, gives 
cos(k y y)e~ kyX ' 



dx' n x (x') 
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(9) 



Before trying to solve this equation, we investigate the long range solution. The nu- 
merator of the dk integral in Eq. (||) approaches zero quickly as x — > oo. Assuming that 
n x (x) ~ as x — > oo, and by setting t = x/x', we must have 1 = --^=r J^ x dt' -k log 



Since dt' j log 
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l-t 
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we can infer that a = - and (3 = 1. This gives 
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n x (x) = — , x — > oo (10) 

for the asymptotic form. 

For the asymptotic form of n y (y), we need to manipulate Eq. (|D|). We first set k' = k\y\ 
which, along with the condition that y — ► oo, enables us to compute the dk integral so that 
we now have n y (y) = — | dx' x , 2 l +y2 . The dx' integral can now be computed by using the 
variable substitution t = x' jy. The final result for the asymptotic form of n y (y) is 

n v(y) = — i — i : \v\ oo (11) 

Eq. (d) can be solved either by successive iterations or by using a matrix method, since 
it is linear. The result is shown in Fig. ||a. n y (y) can be computed from Eq. (|j) since n x (x) 
is now known numerically. The result is shown in Fig. |3|b. 

If we put all the right units back into the result for the n y (y) distribution, we can estimate 
what gate voltage can start the formation of the depletion region on the y axis. This happens 
when the distribution at y = becomes equal to the background density rio y . The condition 
is 

*± = (12) 

In this case, V 9c comes out to be 230 mV. 

This result is not in very good agreement with the experimental resultl. We give two 
reasons to explain this discrepancy. First, our critical voltage value represents the onset of 
depletion, which is not exactly when the source-drain channel is cut off. Electrons can still 
tunnel through if the depletion width is on the order of the Fermi wavelength. Another, less 
important reason, is the behavior of the 2DEG density close to the boundary. The wave 
function for the 2DEG can be represented by ibix.z) = -4f=e kzZ -f- sink x x, where L x and 
L z are chosen boundaries for the gas. Under this assumption, we can calculate the density 
of the gas to be 

n x {x) Ji(2/c j px) 



n 0x k F x 
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(13) 



The effect of this form of charge distribution is to give an effective gate length, which is 
slightly larger than the geometrical one. To get an estimate for the difference, we calculate 
the mean of the extra charge distribution which turns out to be 

Jo u 2 

This is an addition of about 30 A or 0.31 a to the geometrical gate length. We return to 



this again in Sec. IV 



III. DEPENDENCE OF THE DEPLETION ON THE BIAS VOLTAGE 

The Fourier method used in the previous section can not be used when a depletion region 
is present. The Thomas-Fermi equations do not hold in that case. A numerical solution 
had to be employed by discretizing n x (x) and n y (y), converting the integral equations into 
matrix equations, and simultaneously solving for the charge distribution everywhere. 

Similar to the treatment employed in Sec. [H], we first try to investigate the behavior of the 
distribution for large x and y. We assume that n x (x) ~ - and that n y (y) ~ ^ in this limit. 
First, by considering charge neutrality from Eq. (||), we must have J^dx^ + f^dyAr = 0, 
which implies (as I — > oo) that (3 = — ^. Similarly to Sec. |T|, we obtain from the limit 
$x — and 5 y = — ™, which combined with the requirement of Eq. (|7]) gives 

4=4, 6 v = -\ ( 14 ) 

and 

n zO) = ~> x ~> 00 ; n y(y"> = ~2^\y\ ' 00 ■ ^ 

The same results were obtained in Sec. 

To accommodate depletion, only small changes need to be made to Eqs. (|5]) and Eq. (^j). 
Assuming a depletion region from — d/2 to d/2, the new equations are 
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Eqs. ([16]) proved to be a real challenge to solve numerically. The natural logarithms go 
to ±oo when either x or y approaches ±oo or 0. Infinities are also encountered whenever 
x = x' or y = y' . To solve these equations the following technique was employed: A cut-off 
/ was introduced on all the integrals. The integrals now have two main parts. One from 
to / (or — / to /) and another from / to oo (or this and — oo to — /). It turns out that we were 
able to obtain an analytical result to the second portion of these integrals assuming that for 
large enough /, the n x (x) and n y (y) distributions approach their long range approximations. 
The latter parts of these integrals become tails that we attach to the integrals. The attached 
tails are 
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with the latter result achieved by setting t = x/x' in the dx' integral, and t — x/y' in the 
dy' integral. Similarly, we obtain 
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With the addition of these tails and the / cut-offs, Eqs. fllBf ) become 
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The neutrality equation (Eq. (§)) also changes to 
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Eqs. ( |T8"D can now be solved with the additional requirement of Eq. ([7]). A matrix solution 
is obtained by discretizing n x (x) and n y (y) and solving for those values together with S x and 
5 y . The solutions for n x (x) and n y (y) are shown in Fig. |3|. The corresponding solutions for 
Eq. (|7|) are 5 X = 0.749 and 5 y = 0.251, which are in very good agreement with the predictions 
of Eqs. O. 



IV. ADDING WIDTH TO THE ELECTRON GASES 



To get a more realistic result from our simple model, we now try to add width to the 
2DEGs. Assuming that the 2DEGs are in the ground state and extending from —g/2 to 
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^cos(— ). To make the 
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g/2, the wave function for either one is then given by ip{x) = 
numerical calculations easier, we assume a square distribution and pick the width / so that 



f-^dxx 2 I = J^/ 2 dxx 2 ^cos 2 ( 2 ^ 1 ), which gives / = ^Jl — g. The two quantum wells 
for the 2DEGs are 200 A on the x axis and 150 A on the y axis, which give b = 1.268 a® 
and c = 0.951 ao for the x and y axis square width respectively. The new setup is shown in 
Fig. |. 

With the addition of width, Eqs. fllfif) become: 
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The above equations can be solved in a manner similar to that used in solving Eqs. (|T8|). 
The results are shown in Fig. |3|. Fig. ^ shows a plot of the depletion width vs. applied 
voltage for different values of the gate length. Fig. |5] shows a plot of the applied voltage 
needed to create a given depletion width vs. the gate length for various depletion widths. 

We now try to again compare these results with the experimental ones. The experimental 
critical voltage is 450 mV which corresponds to Jr^- ~ 50, shown on Fig. |5]. Assuming an 
effective gate length of about 3 a , this corresponds to a depletion length of approximately 
4.5 \Fy Fig. | shows the voltage distribution on the y axis for various depletion widths. 
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V. AN APPROXIMATE ANALYTIC APPROACH 

Following work done by D. B. Chklovskii et al.i, we make the assumption that the 
potential is constant across the 2DEGs, just like it is in metals. From now on, we will refer 
to this as the metal approximation. The depletion region can be modeled by a continuous 
positive charge distribution. This is shown in Fig. [7], along with all the different boundary 
conditions. To obtain the potential everywhere, we use conformal mapping to transform to 
a region where potential calculations will be easier. In this section, we switch the labels for 
the x and y axis to make the conformal mapping transformation easier. The transformation 
to do this is 

w = sm~ (z/l); z = x + iy, w = u + iv. (21) 

Under this transformation, — I and I go to —tt/2 and 7r/2 respectively, and a goes to a' = 
sinh _1 (a/Z). The boundary condition along the depletion region now becomes 



dv 



<9$ 

v=0+ 



—Icosu. (22) 



D=0" 

We examine the simpler case, where a = 0, which implies that a' in the transformed 
coordinates is now also zero. From symmetry, we need only to consider the region < u < 
Ti / 2. The potential is 

2V 

<$> + (u,v) = —V g + — -u + y A m sm(mu) e~ mv (23) 

7T 



m=evcn 



and 

$_(u,v)= J2 B k cos(ku)e kv (24) 

fc=odd 

above and below the u axis respectively. 

By applying the boundary condition of Eq. (^) and the fact that the two potentials have 
to match at the boundary, we can get equations for the A and B coefficients. The equation 
for the B coefficients is 
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It turns out that the electric field, in the x direction along the boundary, has a singularity 
when x = I (Ref. |9|). This can be canceled by the right choice of I, which is what determines 
the relation between the applied voltage and the depletion length. In the uv coordinates, 
the electric field E x , along the transforms to 



d<5> 

dx 



1 3$ 
= — hm — 

n , u->tt/2 I COS U OU 
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V=0,U-^TT /2 
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V=0,U— >7f/2 



E (-l) ik+1)/2 k (p-) = 0. (27) 



To cancel the 1 / 1 cos u — 1/ y/P — x 2 singularity, we must have — ^ 
gives the condition 

fc^dd v ^ 

A numerical solution is employed to solve for the depletion region width vs. applied voltage. 
The result is 

- = 1%. (28) 

The above equation predicts a linear relationship between the applied voltage and depletion 
length. It also predicts that V 9c = 0. This is shown in Fig. ^. The 2/tt slope is in exact 
agreement with the numerical results. 

Finally, we try to investigate another configuration by using the metal approximation. 

We assume that the transistor has another gate on the other side, so that the new 
configuration is now in the shape of a "cross". This adds more symmetry to the problem, 
which actually makes the solution under the metal approximation easier. We now only have 
to consider the upper half of the uv plane. The new solution for a = a' = is 

v - f- < 29 > 

CIq Hip 

This gives a ir/2 improvement on the slope of the graph of depletion versus voltage over the 
"T" configuration. We also calculate d vs. V g for a ^ 0. These results are shown in Fig. ||. 
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The results are in very good agreement with the high bias limit and predict the correct 
slope between the applied voltage and the depletion width. It should be noted that all the 
slopes converge to the same value in the limit of high bias voltage. Of course, the "cross" 
transistor is not a device that can be manufactured with current technology. However, the 
study of this interesting configuration showed how much improvement an additional mirror 
gate can give. 

VI. CONCLUSION 

The Thomas-Fermi equation allows us to get a good sense of the general behavior of 
depletion formation on the "T" transistor. Because of its simplicity, the metal approximation 
used in Sec. provided a very powerful method for checking our results. Such approximation 
can also be used to get a general idea of a problem without having to explore complicated 
analytical or numerical solutions. And though the method does not predict the correct 
minimum voltage needed to start the formation of the depletion region, it does give the 
correct slope and shows good agreement for large bias voltage. 
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FIGURES 

FIG. 1. Schematic cross section of the T-FET. The two 2DEGs are shown in the gray shaded 
areas. 

FIG. 2. The system of coordinates chosen for the model. We always assume that the positive 
donors are found where the 2DEGs (shaded areas) are. 

FIG. 3. The electron distributions for the gate (x axis) and source-drain (y axis) channels are 
shown for different depletion lengths d and with or without assumed width for the 2DEGs. Lengths 
are in units of ao, and electron distributions are in units of ^^Uq. 

FIG. 4. Plots of the depletion length d as a function of applied voltage V g for different gate 
lengths a. From left to right, the solid lines correspond to gate lengths ^ = 0, 1, 8. The dotted 
line corresponds to the prediction of the metal approximation for ^ = 0. 

FIG. 5. A plot of the applied voltage V g , required to produce a given depletion width d, as a 
function of the gate length a. From bottom to top, the solid lines correspond to depletion widths 

= 0, 5, 10, 15, 20 (5 oo = 1 ^Fy)- The dotted line corresponds to the experimental critical voltage 
of 450 mV. 

FIG. 6. A plot of the potential energy U y (y), in units of the Fermi energy, along the y axis for 
a depletion ^ = 5, 10, 15, 20 (5 ao = 1 A^). The dotted line represents the Fermi level. 

FIG. 7. The boundary conditions for the potential in both the xy and uv coordinate systems. 
The transformation is w = sin~ 1 (z//), with z = x + iy and w = u + iv. The solid lines are the 
positive background charge, and the dotted lines are the 2DEGs, which are assumed to be at a 
constant potential. 

FIG. 8. Depletion length d vs. applied voltage V g for the "cross" transistor using the metal 
approximation for different gate lengths a. From left to right the solid lines correspond to 
-f = 0,1,..., 5. 
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